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Abstract. We give a new proof of Milne's formulas for the number of representations of an 
integer as a sum of 4m^ and 4m(m + 1) squares. The proof is based on explicit evaluation of 
pfafEans with elliptic function entries, and relates Milne's formulas to Schur Q-polynomials 
and to correlation functions for continuous dual Hahn polynomials. We also state a new 
squares. 



1. Introduction 

One of the classical problems of number theory is to count the number of rep- 
resentations of a positive integer n as a sum of k squares. We will denote this 
number by □^(n), where, as is customary, representations 

2 2 
Tl = X-^ ~t~ * ' ' ~t~ X f-. 

that may be obtained from each other by permuting the Xi, or replacing some Xi 
by —Xi, are counted as different. 

The most fundamental results are Gauss' two squares and Jacobi's four and 
eight squares formulas: 

d\n, d odd 

□^(n) = 8 ^ ci, (1.1b) 

d\n, i\d 

□sH = 16 Y^{-iY+^£, (l.lc) 

d\n 
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where the sums run over positive divisors. These may be compared with Legendre's 
formulas for sums of triangles: 

d|4n+l 

A4(n) = XI ^' ^^-2^) 

d\2n+l 

d\n+l, {n+l)/d odd 

It is known that, in general, ^2k{n) can be written as the sum of two terms, the 
first being an elementary divisor sum and the second the n-th Fourier coefficient of 
a cusp form. Moreover, the second term vanishes only for /c < 4, which indicates 
that there is no very simple extension of (1.1) to more than 8 squares. 

A novel approach, motivated by ajfine superalgebras, was initiated by Kac and 
Wakimoto [KW] , who used denominator formulas for such algebras to derive sev- 
eral new infinite families of identities, both for squares and triangles. A particu- 
larly interesting case is the "queer" series of affine superalgebras Q{m), for which 
the denominator formula was merely conjectured. For Q(2m — 1) and Q{2m), 
respectively, it implies the triangular number identities 

^'"'■(") = 4„(.„-.)n'"-,. E n «-*|)^ (1.3a) 

^^j=^ Aii/iH ^kmlm=2n+m^ ^=^ l<i<j<m 

ki>k2>--->km 
ki and li odd positive 

A4^(„+i)(n) = -^ n (^'-^D'- (l-3b) 

^^j='^-^' klh + ■■■+k,-„lr,^=n+\rn(:m+l) *=1 l<i<j<m 

kl>k2>--->krn 

ki positive, li odd positive 

The case m = 1 gives Legendre's formulas for four and eight triangles. The sums 
of squares formulas contained in [KW] are different from those discussed in the 

present paper. 

The identities (1.3) were first proved by Milne [M3], see also the research an- 
nouncements [Ml, M2], using an approach different from Kac and Wakimoto. 
Independently, the more general denominator formula was proved by Zagier [Z]. 
Getz and Mahlburg [GM] showed that it contains further triangular number iden- 
tities, such as 

A2m{n) = Yl (-l)lft'=^-i-2^ (1.4) 

fei/iH \-kmlm=4:rnn+ivP 

ki and li odd positive 
kiS±{2i—l) mod 4m 

which reduces to (1.2 a) for m = 1. 
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In [Rl], we made the observation that the denominator formulas for queer affine 
superalgebras can be written as pfaffian evaluations, for matrices with eUiptic func- 
tion entries. Moreover, we showed that they follow from a classical determinant 
evaluation due to Frobenius [F]. We also attempted a complete study of the 
implied triangular number identities, finding formulas for Am? /d triangles, when 
d I 2m, and Am{m + 1)/ d triangles, when d \ 2m or d \ 2m + 2. As an example, let- 
ting 0? = 2 in [Rl, Eq. (4.4b)] gives a 2m^ triangles identity, which can be written 
as 

lij=l \J-> kih+--+k„,l^=An+'m? 

kl>k2>--->km 
ki and li odd positive 



n {{-i)'^^'^-'^k,-{-iy^^'^-'^k?j\ (1.5) 



X 

l<i<j<m 

Milne's proof of (1.3) also uses elliptic functions, but is based on continued frac- 
tions and Hankel determinants, rather than on pfaffians. With similar methods 
he obtained new formulas for sums of 4m^ and 4m(m -|- 1) squares, see Corol- 
lary 4.5 below. As an example, Milne's 16 squares formula can be written [M3, 
Corollary 8.1] 

kl=n, k,l>l 

+ ^ J2 {-if^^-'^^'^-'^+^'^^-'^^'^-'^kMkl - klf. 
o 

kih+k2l2=n 

fel>fe2>l,h,/2>l 

For a readable introduction to Milne's work, we refer to the survey [CK] . 

Ono [On] derived seemingly different formulas for 4m^ and 4m (m + 1) squares 
from (1.3), using the fact that the generating functions for squares and triangles are 
related by a modular transformation, cf. (2.9) below. Moreover, in the introduction 
he indicates an alternative proof, also using modular forms but without relying on 
(1.3). In the Appendix, we will show that Ono's formulas are in fact equivalent to 
Milne's, although his proof is different. Another modular forms proof of Milne's 
formulas was recently given by Long and Yang [LY] . This should be related to the 
alternative proof indicated by Ono. 

The purpose of the present work is to extend the analysis of [Rl] from triangles 
to squares, deriving Milne's sums of squares formulas from elliptic pfaffian evalu- 
ations. These evaluations follow from Frobenius' classical determinant, and may 
be viewed as modular duals of those used in [Rl]. Thus, just as in Ono's proof 
square number identities arise as modular duals of triangular number identities, 
although the details are different. 
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To be precise, wc give two, closely related, derivations of Milne's formulas. Ap- 
plying two different modified Laurent series expansions to the matrix elements 
of our pfaffians, we obtain two different analogues of the Kac-Wakimoto denom- 
inator formula. Theorem 3.5 and Theorem 5.5. Theorem 3.5 looks very similar 
to the Kac-Wakimoto formula, though convergent Lambert series are replaced 
by Abel means of divergent series. Theorem 5.5 involves convergent series, but 
Schur polynomials are replaced by Schur Q -polynomials. To be precise, classical 
Schur Q-polynomials are labelled by positive integer partitions; here, we need an 
extension to the case when some labels are negative. 

Theorems 3.5 and 5.5 involve a number of free variables. Speciahzing all vari- 
ables to 1 gives formulas for Am? and Am{m -\- 1) squares. The resulting formu- 
las are equivalent, though to see that is far from obvious. Starting with Theo- 
rem 3.5, one readily obtains the Hankel determinant form of Milne's identities. 
Corollary 3.6. Applying an identity relating Hankel determinants to correlation 
functions of orthogonal polynomial ensembles, we deduce explicit sums of squares 
formulas involving correlation functions for continuous dual Hahn polynomials, see 
Corollary 4.4. Expanding the correlation functions into Schur polynomials yields 
the Schur function form of Milne's identities. Corollary 4.5. 

On the other hand, specializing all variables in Theorem 5.5 to 1 yields sums 
of squares formulas involving similarly specialized Schur Q-polynomials, Corol- 
lary 5.7. Such quantities, and more generally Schur Q-polynomials specialized to 
a geometric progression, are studied in [R2]. It follows from that work that the 
Schur Q-polynomials appearing in Corollary 5.7 agree with the correlation func- 
tions of Corollary 4.4, and consequently that Corollary 5.7 is equivalent to Milne's 
identities. 

It is natural to try to generalize Milne's formulas to other numbers of squares, 
similarly as was done for triangular numbers in [CM, Rl]. In principle, it should 
be possible to obtain formulas for Aw? /d squares, when d \ 2m and Am{m + l)/d 
squares, when d \ 2m or d | 2{m + 1), from our elliptic pfaffian evaluations. 
However, in practice the computations are difficult to handle, and the resulting 
identities seem very complicated to state. As an indication of what the results 
may look like, we state without proof a new formula for 2m^ squares in Theorem 
6.1. In particular, this embeds the two squares formula (1.1a) in an infinite family, 
similarly as Milne's identities do for four and eight squares. 

The reader may question the merits of our new proof of Milne's formulas, com- 
pared to the simple modular forms proofs of Ono and of Long and Yang. First 
of all, in our approach (similarly as in Milne's original proof) the identities are 
derived from classical results. Although modular forms are a powerful tool for ver- 
ifying Milne's formulas, they seem less useful for deriving the results from scratch. 
Second, our proof explains the coefficients in Milne's formulas (i.e. the double 
sums over A and /j, in Corollary 4.5), by relating them to Schur Q-polynomials and 
to correlation functions for continuous dual Hahn polynomials. Third, we obtain 
Milne's formulas as a special case of the more general Theorems 3.5 and 5.5, which 
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might be of independent interest, plausibly in the theory of super algebras. In con- 
clusion, we feel that the three known approaches to Milne's identities (Milne's 
original proof using continued fractions and Hankel determinants, the modular 
forms proofs of Ono and of Long and Yang and the present work using pfaffian 
evaluations) each add something to the understanding of these deep and beautiful 
results. 

Acknowledgements: I thank Stephen Milne for several useful comments. 

2. Preliminaries 

2.1. PfafHans. As in [Rl], we define the pfaffian of a skew-symmetric matrix as 

1 



M 



^pfaff^(a,) = 



'^cr(2i-l),(7(2i)) 



where M is the integral part oi 111/2. This definition is standard when m is even, 
but not when m is odd. In even dimension, 

pfaff (ttij) 

K.l<i,j<2m 



i.det {aij) 

l<i,3<2m 



(2.1) 



The cases of odd and even dimension are related by 



pfaff {ttij) = pfaff 

l<i,j<2m+l l<jJ<2m+2 



V-i 



-1 0/ 



We note in passing the pfaffian evaluation 



pfaff 

l<i,j<m 

\Xi + X 



= n 



X j X A 



l<i<j<m + 



(2.2) 



When m is even, this is a classical identity of Schur [Sc]. The case of odd m then 
follows by letting x^ — O. 

2.2. Theta functions. Throughout, q will be a fixed number with < g < 1. 
We will use the notation 



(a)oo = (a; q)oo = JJ(1 - aq''), 



k=0 



(ai, . . . , ttrnjoo — (fll, • • • 5 Cim', q)oo — {<ll'i Q)oo " " " {(^m', q)oo- 

When the base is suppressed from the notation, it is always taken to equal our 
fixed number q. 

We introduce the theta function 



9{x) = 9{x; q) = {x, q/x; q)oo- 
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It satisfies 

e{x-^) = e{qx) = -x-^e{x) (2.3) 

and the modular transformation 

Tlie Laurent expansion of 6 is given by Jacobi's triple product identity 

00 

{qUe{x)^ J2 (-1)^®^^ (2.5) 
fe=— 00 

We also mention the Laurent expansion 

iq)l,0iax) ^ x'' 



Yl 1 fe' 5 < kl < 1, (2.6) 

1 — aa'^ 



0{a)9{x) 1 - aq 

which is a special case of Ramanujan's iV'i sum [GR, Eq. (5.2.1)], together with 
its limit case 

9'ix) x^ , , . „x 

-^ = -Et3^. 9<W<1- (2-7) 

2.3. Generating functions. The triple product identity (2.5) implies explicit 
formulas for the generating functions for squares and triangles. Indeed, let 

00 00 



□(?)= E ?"' = l + 2$]g"', 

n=— 00 n=l 

2 

n=— 00 n=0 

so that, according to the standard conventions that we use, 

00 

□(g)^ = J]nfe(n)g", 



n=0 

00 



n=0 

Then, by (2.5), 

□(g) = (g'; g')oo e{-q- q^) = {q\ -q, -q; q^). 



9 9\ 



A(g) = ^ (g; g)oo ^(-g; q) = {q, -q, -q; q) ^'^ ' ^ 



2 (g;g 



2^ 

00 
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Note that the special case x = 1/2 + i/h oi (2.4) gives 

A(e-2^/'*) = ^ e""'^^ □(-e-^'^). (2.9) 
We recall the Lambert series versions of (1.1), that is, 

oo ^ 

fe=l ^ 

fc=i ^ 



□(#=1 + 16^^^^. (2.10c) 

Expanding the denominators as geometric series leads immediately to (1.1a) and 
(1.1c). In the case of four squares, one obtains 

8j](_i)(<^-i)(i-i)d, 

d\n 

which is equivalent to (1.1b); cf. the case /c = 1 of (A.l) below. 

In §6, we need two further identities of Jacobi. First, we have the Lambert 
series expansion 



,Atf)*n(,y^ = $:^. (2.11) 



Second, comparing (1.1b) and (1.2b) gives 

8A4(n) = 04(271 + 1), 
which can equivalently be written 

16gA(gY = □(#-□(-#. (2.12) 

2.4. Sums of triangular numbers. For comparison, we briefly recall how the 
triangular number identities (1.3) follow from the pfaffian evaluations 

^i^^ (-^pI^\^q\rnirn-l)YX^rn-i TT V^^) (2.13a) 
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For further details we refer to [Rl], though the essential points are contained 
already in [Z]. 

It will be convenient to introduce the functions 

deti<i i<m(xt') 

S,{x,,...,xm) = ^ i<M<M^; ^^^^^ ^2.14) 

which are essentially Schur polynomials. Indeed, since is anti-symmetric in the 
variables /li we may assume ni > ■ ■ ■ > Hm- Moreover, since 

S{i^i+a,...,ij.m+a){^lj • • • ) ^m) = ^1 ' ' ' S{ni,...,Hm){^lj • • • ' ^m) (2-15) 

we may also assume /i^ > 0. Then, writing /i^ — + m — i, so that Ai > • • • > 
Am > 0, Sfj, equals the Schur polynomial s\. 

Returning to the task at hand, we expand each matrix element in (2.13) as a 
Laurent series in the annulus ^yq < \xj/xi\ < Applying (2.6) and (2.7) one 

arrives after some elementary manipulation at the multivariable Lambert series 



(g)^ TT 1 TT iqXj/xi,qXi/xj)^ _ 



1=1 » l<i<j<2m 

m 1 1.. 



IT:; 7— T 'S'(A;i,...,fc^ _fci_i)(a;i, . . . ,a;2m), (2.16a) 

1 ^„ 1 1 — O '"^2 



fel>-->fem>0 i=l 



(g)^ TT — TT {qXj/x,,,qx,-Jxj)oo 

(V^)^ ii i<i<}<2m+l (V^^i/^i' \/^Xi/Xj)oc 

m If,. 

^ YlYir^,S^'^^,-M-k^^^^^^^ (2.16b) 

A;i>->A;ru>l i=l 

which are equivalent to the denominator formula for Q{2m — 1) and Q{2m), re- 
spectively. 

Specializing Xi = 1, the left-hand sides of (2.16) reduce to A(y^)^"^ and 
^(V^)^"^^'"^^\ respectively. On the right, applying the classical formula 

S,{n= n (2.17) 

and expanding the denominators as geometric series leads after simplifications to 
(1.3). 

More generally, one may let Xi — lo'^ ^, with u a suitable root of unity. This 
leads to more general triangular number identities such as (1.4) and (1.5), see [Rl]. 
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2.5. Hankel determinants. We recall the following classical result [I, Corollary 
2.1.3]. 

Lemma 2.1. Let 

f^if) = / f{x)dn{x) 



he a linear functional defined on polynomials, let Ck — ^jl{x^) be its moments, and 
let 

A(x)= W {x^-Xi)^ det ix>-^) 

^ ' -L-L ^ ' l<i,j<m ^ * ^ 

l<i<j<m 

denote the Vandermonde determinant. Then, 

det {ci+j-2) = I ^(a^i, • • • , Xmf dii{xi) ■ ■ ■ dii{xm)- 

l<t,j<m ml J 

Proof. Clearly, 

/m 
A{xi, Xm)^ dll{xi) ■ ■ ■ dll{Xm) = ^ Sgn(c7) Sgn(r) C^(i^+r{i)-2- 

(T,TeSm i=l 

Replacing a hy ar, this may indeed be written 

m 

Y] sgn(a) TT Ci+a{i)-2 = m! det (q+j-s). 



(j,Tes, 



□ 



Equivalently, defining the coefficients C{ki, . . . , km) by 

m 

n i^j-^if= X] C{ki,...,km)Ylx'l\ 

l<i<j<m k\,...,km *=1 

we have, for arbitrary scalars c^, 

m 

V C(A;i, ...,km)T\c),^^m\ det (^+^-2) . 

This version of Lemma 2.1 will be useful for our discussion of Ono's identities in 
the Appendix. 

If pk and qk are arbitrary monic polynomials of degree k then, by linearity, 

^4et^(/"fe-i?j-i)) = ^Jet^ [l^{x'+''^)) . 

In particular, if /i is a positive functional we may choose pk = as the corre- 
sponding monic orthogonal polynomials. Then, Lemma 2.1 gives 

— / A(xi,...,xj2d/i(xi)---d/x(xJ =TT||pi_i||2. (2.18) 

i=l 
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2.6. Tangent numbers. We will need the following evaluation of Abel means of 
alternating power sums in terms of tangent numbers, see [T, Theorem 2.5]. It is 
equivalent to the classical evaluation of Riemann's zeta function at the negative 
integers, though for completeness we provide a self-contained proof. We find it 
convenient to define tk by 

tan — = > ti;— -T. 

2 ^ (2A;-1)! 

Equivalently, in standard notation for tangent and Bernoulli numbers, 

, ^i-2kr (4'-l)l^2fe| 

tk — Z 1^ — — 



k 



We will denote Abel means of possibly divergent series by 

Ck = lim ^t'^'cfc, A C Z. 



keA 



Lemma 2.2. One has 

oo 

E'(-i) 



k=l 



1/4, m = 0, 

0, m = 2,4,6, . 

(-l)"+ii„/2, m = 2n - 1. 



Proof. We assume that m > 0. Since 

1-x 



we may write 



2^(-l)'=+^A;"^ 



1 + x 



l + 2^(-l)V, |a;|<l. 



k=l 



k=l 



X 



dx 



x=l 



1 — ,T 
1 + X 



t=0 



1 - e* 
1 + e* 



t=o 



i tan ■ 



2i 



which implies the desired result. 

We introduce the moment functionals 



0, 1. 



k=l 



CoroUciry 2.3. The moments of jiQ and /ii are given by 



k-U 



(2.19) 



(2.20) 



(2.21) 



□ 



(2.22) 



(2.23) 



l2o[X ) = j^i[X 

In Lemma 4.1, we give alternative expressions for the functionals /i^. As a conse- 
quence, we shall see that the corresponding orthogonal polynomials are continuous 
dual Hahn polynomials. 
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3. From pfaffians to Hankel determinants 

We will obtain sums of squares formulas from the following pfafiian evaluations, 
the first of which was already given in [Rl]. 



Lemma 3.1. One has 



pfaff 



9{xj/xi) 



n 



d{Xj/Xi) 



l<i7<2m \0{-Xj/Xi) J ^<.^^:^2^ 0{-Xj/Xi) ' 



(3.1a) 



pfaff I 1 

i<i,j<2m+l 



^ Xj 9 ( Xj/Xj) 
Xi Xj/ Xi^ 



(9) 



2m 

oo 



n 



9{xj/xi) 



(3.1b) 



Remark 3.2. Since 6{x; 0) — 1 — x, Lemma 3.1 reduces to (2.2) when g = 0. A 
different elliptic extension of Schur's pfaffian evaluation was recently obtained by 
Okada [O], see also [Rl, Remark 2.1]. 

We give two proofs of (3.1b). The identity (3.1a) can be proved similarly [Rl]. 
First proof of (3.1b). We start from the Frobenius-Stickelberger determinant [FS] 



det 

l<i,j<n+l 



xjyi O'jxjyi) 

V-1 ••• 



1\ 

1 



• -1 0/ 

_ g(xi • --Xnyi ■■■vn) Ui<i<j<n^jyA^i/^jWiyi/yj) 

Choosing n — 2m + 1 and yj — —1/xj, and moreover adding 1/2 (gf)^ times the 
last column to all other columns gives 



det 

l<i,j<2m+2 



1 



-1 



1 + 2 



Xi Xj I Xj^ 



1 

0/ 



n 



9{xi/xiy 



We now observe that the matrix on the left is skew-symmetric. This is easily 
proved by differentiating (2.3). Thus, in view of (2.1), we may conclude that 
(3.1b) holds up to a factor ±1, which is independent of q by continuity. By 
Remark 3.2, that factor has to be +1. □ 
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Second proof of (3.1b). In (2.13b), let q = e and 
logarithmic derivative of (2.4), that is. 



Applying the 



2ie 



2_^'(e2™;e-2-A) 



0{e 



2-Kix- , 



-2% 



/h) 



i + h — 2hx — 2he 



-2-Khx 



-2-Khx. g— 27r/i^ 
2-Khx. g-27r/i^ ' 



the left-hand side of (2.13b) takes the form 

'hp^/^\^ / 2-Khzj nt( _ 2-Kh{zj-Zi). -2-nh\ 

2i 



+ 2{zj - Zi) 



By Lemma 3.3 below, we may subtract 2[zj — Zi) from each matrix element, which 
gives a pfaffian of the desired type. Similarly applying (2.4) and (2.9) to the 
right-hand side of (2.13b) completes the proof. □ 

The following lemma was used above. 

Lemma 3.3. For any odd- dimensional skew- symmetric matrix {aij), 

pfaff {aij -\- bi — bj) = pfaff (%). 

l<i,j<2m+l l<i,j<2m+l 

Proof. Adding multiples of the last row and column to the previous ones gives 



det 

l<i,j<2m+2 



1\ 



dij + bi- bj 



V-1 



1 

-1 0/ 



det 

l<i,j<2m+2 



( 



V-1 



"13 



1\ 

1 

-1 0/ 



By (2.1), we may extract square roots to conclude that 

pfaff (ojj -|- 6i — 6j) = ± pfaff (oij). 

Since both sides are polynomial expressions, the sign may be determined by letting 
6i = 0. □ 

Proceeding in analogy with §2.4, the next step would be to expand the matrix 
elements in (3.1) as Laurent series in the variables Xi. An important difference 
from (2.13) is the presence of singularities at Xi — —Xj, which precludes Laurent 
expansion near Xi = Xj = 1. To circumvent this difficulty we will use Abel means, 
for which we recall the notation (2.20). 

Lemma 3.4. // \x\ — 1 and x ^1, then 



9' 



1 + 2x 



X] 



e{-x) 



(3.2a) 
(3.2b) 
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Proof. By (2.6), 

j^mi. ,<i.i<, ,3,3, 

Replacing x by qx and using (2.3) gives 

WMf) . f , (3.4) 

Letting f{x) denote the left-hand side of (3.2a), we write 

t — >^ Zi 

Expanding j{xt) using (3.3) and f{x/t) using (3.4) yields (3.2a). 
As for (3.2b), if f{x) denote its left-hand side then (2.7) gives 

' x)^ 



It is easy to check that f{l/x) = —f{x). Indeed, this is equivalent to the skew- 
symmetry of the matrix (3.1b). Thus, 



Similarly as above, we can now deduce (3.2b). □ 
Applying (3.2a) to the matrix elements of (3.1a) gives 



/(.) = -l + 2^^^^, i<N<,-. 



l<i<j<2m ^ ■>' 



kl,...,km¥^0 i=l ^ aeS2m «=i ^ > 

where we assume that \xi\ — \xj\ and Xi ^ —Xj for all i and j. Multiplying both 
sides with nj<j(^i + ~ ^j) adopting the notation (2.14), this may be 
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written 



l<i<j<2m \ ^ J' " ^ J/^ l<i<j<2m 



- + g 

Since the summand is anti-symmetric and even as a function of ki, and vanishes 
unless the ki arc all distinct, we may restrict the summation to ki >■■■> km ^ ^ 
if we multiply by 2™m!. Moreover, 

This gives the first half of following result, which should be compared both with 
the Kac-Wakimoto identities (2.16) and with the closely related expansion given 
in Theorem 5.5 below. The second half follows similarly from (3.1b). 

Theorem 3.5. // \xi\ = \xj\ and Xi ^ —Xj for all i and j, then 



(qT^ tt {QXj/Xj, qXj/Xj)^ ^ -pr . X 



m ^ u. 



fcl>-->fcm>l 1 = 1 

(-?)^ l<j<j<2m+l -QXi/Xj)oo i<i<j<2m+l 



To obtain sums of squares formulas from (3.5), we let Xi = 1. Then, the left- 
hand sides reduce to □(— g)^™ and □(— g)'^™('"+i). On the right, (2.17) gives 

S(ku-,km-km,--ki) (1^™) = ^2m-l ., II n ^^i ~ 

lii=l ^- i=l l<i<j<m 

'S'(fci,...,fc^,0 -fern,... -fel)(l^"''^^) = p.2m ., II ^^i II (^i ~ ^If- 

lli=l ^- i=l l<i<j<m 

Focusing on (3.5a), we thus obtain 

- y TT(_i)..+il-(-g) 
I 11^ ^ l + (-q)ki 

'ki,...,km=l j=l ^ ^ l<i<j<m 



"*-lii=l fc™=li=l -"^^ '''' l<i<-i<m 
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By Lemma 2.1, this may be written in Hankel determinant form as 

om(2m— 1) 



where 



Next, writing 



leads to the decomposition I'o — Xo + Ho, where /Xq is as in (2.22), and 



since the latter sum is convergent, there is no need to take the Abel mean. Us- 
ing Corollary 2.3 to identify the moments of /iq, we arrive at the first half of 
Corollary 3.6. 

Similarly, (3.5b) leads to the identity 

om(2m+l) 

where 

1 + (-?)' 



Writing i^i — Xi + fXi, with fj,i as in (2.22) and 

k=l ^ ^' 

yields the second half of Corollary 3.6. These identities are equivalent to [M3, 
Theorems 5.3 and 5.5]. 

CoroUairy 3.6 (Milne). One has 

9m(2m-l) / °° „feL2j+2jr-3' 

□(g)4^^ ^ . , det \u+., + A{-iy+^y ^\ , 

9m(2m+l) / °° fcr,2j+2j-l~ 

□ ( )4m(m+l)^ ^ det t,+,+4(-l)^+^ y — 

where the numbers tk are defined in (2.19). 

Let Pk\x', q) be the monic orthogonal polynomials corresponding to u^. Then, 
(2.18) gives the following reformulation of Milne's formulas. 
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Corollary 3.7. In the notation above, 

9'm(2m-l) ^ 

11^=1 ^! f=\ 

Om{2m+l) 

Equivalently, 

\\pf\^;,)r='^^i^n(,r-\ (3.6a) 

llpi"(^;.)f^^?^if^ (3.6b) 

It is well-known that any positive definite Hankel determinant can be expressed 
in terms of orthogonal polynomials. The point here is the explicit expressions 
for the moment functionals, which might suggest that the polynomials p^^\x;q) 
are of independent interest. In particular, an alternative proof of (3.6) would 
lead to a new proof of Milne's sums of squares formulas. Note also that p^^\x; q) 
are g-analogucs of the polynomials p^i^\x;0), which, as we shall see in §4.1, are 
continuous dual Hahn polynomials. However, they are of a different type from the 
polynomials in the g-Askey Scheme [KS]. 



4. Correlation functions and sums of squares 

To obtain sums of squares formulas from Corollary 3.6, we must expand the 
determinants as power scries in q. In [M3], this is achieved using the Cauchy-Binet 
formula, which leads to identities involving Schur functions, see Corollary 4.5. We 
will give a more conceptual derivation of these identities by relating them to 
correlation functions for continuous dual Hahn polynomials. 

4.1. Continuous dual Hahn polynomials. We need the following alternative 
expressions for the functionals defined in (2.22). 

Lemma 4.1. One has 

coo _2£+l X/'™2^ 



Proof. Write 



1 /'°° r^^+Vfr^l 



°° x^'+^f{x^) _ f°° x^'+^e'^'f{x^) 



Jo SinhfTTX) A-^0+ J_ 



2 7q sinh(7rx) a-^o+ sinh(7ra;) 

and expand as the sum of residues in the upper half-plane. □ 
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Alternatively, Lemma 4.1 can be deduced from Corollary 2.3. One is then 
reduced to the integral evaluation 

^ x^''-^ . t 



Iq sinh(7ra;) 2 
which can be found in standard tables, or derived from 



^0 



sinh(tx) ^ tan(V2) _^ ^ , ^ ^ 
smh(7rxj 2 



by expanding both sides as power series in t. 

Let (pjf^)^o monic orthogonal polynomials corresponding to ji^, s — 0,1. 

We claim that they may be identified with continuous dual Hahn polynomials [KS] . 
In general, when a,b,c> 0, one denotes by {—l)''Sh{x; a, b, c) the monic orthogonal 
polynomials with respect to 



/ 



r{a + ix)r{b + ix)r{c + ix) ^ 



r{2ix) 



fix"^) dx. 



Using elementary properties of the gamma function, it follows from Lemma 4.1 
that 

pt\x) = i-lfSkix; 0, 1/2, 1), p^^\x) = {-lfS,{x; 1/2, 1, 1). 
Since p^^\x) — Pk\x', 0), the case = of Corollary 3.7 reads 

m 2m— 1 m ^ 2m 

1=1 i=l i=l 1=1 

which agrees with the known expressions for the norms [KS]. 

We remark that and p^j^^ combine naturally to a single orthogonal system. 
Indeed, if 

P2k{x) = Pk\x'^), P2k+i{x) = xp'j^\x'^), 
then ipk{x))'kLo are the monic orthogonal polynomials corresponding to 

xf{x) 
) sinh(7rxj 

These are Meixner-PoUaczek polynomials; in the notation of [KS] 
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4.2. Correlation functions. We recall some general facts about correlation func- 
tions, see [R3] for references and further details. 
Let 



Kf) = I fix) di^ix) 



be a positive moment functional and (pfe(a;))^o the corresponding family of monic 
orthogonal polynomials. Normalizing 

A(a;i, . . . , dfi{xi) ■ ■ ■ dfj,{xm) 

to a probability measure, it defines an orthogonal polynomial ensemble. Such 
ensembles arise in several contexts, including the theory of random hermitian 
matrices [K] . An important role is played by the correlation functions 



/ 



A{xi, Xnf dl^{Xm+l) ■ ■ ■ dl^{Xn), < Ul < U. 

We will normalize these functions as 

C[X) — C [Xi, . . . , Xjn) — 7 M T-rn m mo a / \2 

X J A{xi,...,XnfdlJ,{Xm+l)---dlJ,{Xn). (4.2) 

Then, C is a symmetric polynomial. 

There are several expressions for correlation functions in terms of orthogonal 
polynomials, including 

(^/^N _ 1 / Pn{Xi)Pn-liXj) - Pn-l{Xi)Pn{Xj) \ 

\\Pn-l\\'^"'A{x)^ l<i,j<m\ Xi-Xj J ^' ^ 

- ^det 1^^^-' I (4.3b) 



Kl ||Pn-^||2A(x)4 l<^,j<2m ^ [p'^_^^j_,{Xi) , TU + I < I < 2m ^ 

_ 1 (deti<ij<m(pfc,(a;j))) 



V- (deti<ij<m(Pfc.(a^i)))^ 
where the diagonal entries in (4.3a) are interpreted as the limit 

,. Pn{Xi)Pn—l\Xj ) ~ pn—l{Xi)Pn{Xj ) I / \ / \ I / \ / \ 

lim ^ Pn{Xi)Pn-l{Xi) -Pn_i{Xi)Pn{Xi). 

To recover Milne's Schur function expansions we need another identity, namely 
[R3, Proposition 1.8], 

^ (_l)EI^i(Ai+Mi) 

n^x ^ lli=l ^e[n]\S,3e[n]\T 

0<Xm<- -<M<n—m ^-^i-—^ " 

0</Um<---<;ii<n— m 
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where Ck = and where 

S — {Afe + m + 1 — k; 1 < k < m}, T = {/i^ + m + 1 — k; 1 < k < m}. 

An important relation between Hankel determinants and correlation functions 
follows from Lemma 2.1 upon replacing fihj ^ + X, where we think of /i as positive 
and A as arbitrary. Since the integrand is symmetric, we may write 

det U(x'+^-'') + A(x'+^-')) 

\ / A(a;i,...,x„i)^ci(/x + A)(a;i)---(i(// + A)(x^) 
m! J 

""if 

y] -77 TT / A(xi, . . . , Xmf dX{xi) ■ ■ ■ dX{xs)diJ,{xs+i) ■ ■ ■ dfJ,{Xm), 

s\{m-sy. J 

where the integral over the last m — s variables is a correlation function. We thus 
arrive at the following result, which is a standard tool of random matrix theory; 
see e.g. [J, §2]. 

Lemma 4.2. In the notation above, 
det + A(x'+^-')) 

l<i,i<m ' 
m m „ 

= nil^^-ill'El / A(xi,...,x,)2cr(xi,...,x,)dA(xi)-.-dA(x,). 

4.3. Milne's sums of squares formulas. By Lemma 4.2, the Hankel deter- 
minants of Corollary 3.6 may be expressed in terms of correlation functions for 
continuous dual Hahn polynomials. If we let C^'^ denote the correlation function 
defined by choosing = /i^ in (4.2), then applying Lemma 4.2 with /i — /iq and 
A = Ao yields, using also (4.1), 



m 



DW""" = E7 E Utti^ n «-*?)^CT°(-*?,...,-a. 

s=0 ki,...,ks=l i=l ^ l<i<j<s 



It is understood that the sum over ki, . . . ,ks equals 1 if s = 0. Similarly as in 
Theorem 3.5, we may restrict the summation to /ci > • • • > /c^ if we multiply by 
s!, thus obtaining the first half of Corollary 4.3. The second half follows similarly, 
choosing /i — /ii and A = Ai in Lemma 4.2. 

Corollary 4.3. One has 

m s k J 

nfer"= = E*' E Uttt^ n «-*|)=cT»(-*;,...,-*s, 

s=0 ki>- >ks>l i=l ^ ^' l<i<j<s 
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□ (^g)4m(m+l) 



s=0 ki>->ks>l i=l ^ ^' l<i<j<s 

Expanding 

^ oo ^ oo 

gives the following sums of squares formulas. We have excluded the terms with 
s = 0, since they only contribute to the trivial coefficients corresponding to n = 0. 

Theorem 4.4. For n > 0, 



kili+-- 

ki>-: 
h,.-, 

X C^'^{—kl, . . . , —kg), 



s=l kili-\ \-ksls=n i=l l<i<j<s 

ki>->ks>l 

h,...,ls>l 



s=l fei^iH \'kglg=n i=l l<i<j<s 

ki>--->ks>l 
h,...,ls>l 

X C^'^{—kl, . . . , —kl). 

Expressing the correlation functions as in (4.3) gives explicit versions of these 
sums of squares formulas in terms of continuous dual Hahn polynomials. Using in- 
stead (4.4), together with the expressions (2.23) for the moments and (4.1) for the 
norms, yields the following identities, which are equivalent to [M3, Theorems 7.1 
and 7.2]. 

Corollary 4.5 (Milne). For n > 0, 

nm(2m-l) ™ « 



=1 kihA \-ksls=n 1=1 l<i<j<s 

ki>--->ks>l 

h,...,ls>l 



^ E . , , det (ti+j_i)sA(A;^,...,A;2)s^(/c^,...,/c^), 

0<Xs<-<Xi<m—s 
0<iJ,s<-<lJ-i<m—s 
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9m(2m+l) « 
lii=l^- s=l |_j,^;^=„j=i l<i<j<s 

fcl>--->fcs>l 
h,...,ls>l 

0</Us<---</ii<m— s 

where is as in (2.19) and 

S ^ {Xk + s + 1 - k; 1 < k < s}, T ^ {^k + s + 1 - k; 1 < k < s}. 



5. Relation to Schur Q- polynomials 

In [R2], we showed that the correlation functions appearing in Theorem 4.4 also 
arise from generalized Schur Q-polynomials by specializing all variables to 1. In 
this section, we give a slightly different proof of Milne's formulas, where Schur 
Q-polynomials appear naturally. 

5.1. Schur Q-polynomials. When m < n and A e Z'", we write 



Qx{xi,..., 




l<i<j<n 



Here, S'„ acts by permuting the variables Xi,....Xn and Sn-m is the subgroup 
acting on Xm+i, ■ ■ ■ ,Xn- When Ai > A2 > • • • > > 0, this is a standard 
definition of Schur Q-polynomials [Ma]. Note that there is no analogue of (2.15), 
so the case when some Aj are negative cannot immediately be reduced to the 
classical situation. 

We will not work with (5.1) directly, but rather with the alternative expression 
[R2, Lemma A.l] 

Q\{xi, ...,Xn) 

l<i<i<n aGS„ i=l i=l •^<y(m+2t-l) ^<T(m+2i) 

where k is the integral part of (n — m)/2. In analogy with (2.14), the right-hand 
side can be rewritten as a quotient of two pfaffians [N, Eq. (A 12)]; however, the 
form given here is more convenient for our purposes. 

One of the main results of [R2] is the identification of Qx{l,q, ■ ■ ■ with 
generalized Christoffcl-Darboux kernels for continuous g-Jacobi polynomials. In 
the case when q = 1 and A = (Ai, . . . , A^, — A^, • • • , — Ai), the continuous g-Jacobi 
polynomials reduce to continuous dual Hahn polynomials and the generalized 
Christoffel-Darboux kernels to correlation functions. To be precise, noting that 
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the kernels K^^"^ featuring in [R2] are related to the kernels C^'^ introduced in §4.3 
through 

it follows from [R2, Corollary 5.11] that 

2n\ 



Q(Al,.--,Am,-Amr--,-Al)(l ) 



= 8-n^^ n i>^^->^]fC:f{-Xl...,-Xll (5.3a) 

i=l l<i<j<m 

<5(Ai,...,A,„ -A^,... -Ai)(l^"^^) 

m 

=(-i)-8-n^' n (5-3b) 

i=l l<i<j<m 

5.2. A second proof of Milne's formulas. Wc return to the pfafRan evalua- 
tions (3.1). In §3, we handled the singularities at Xj = —Xj by introducing Abel 
means. We will now use a slightly different approach, removing the singularity 
by subtracting a rational function. This corresponds precisely to the previous 
decomposition /i^ — + of functionals, so the two ideas are in fact intimately 
related. 



Lemma 5.1. For q <\x\ < q ^ , 



1 + 2x^ = + 2 f i-'in^' - , (5.4b) 

e(-x) l + x ^ l-q^ ^ ' 

fc=i 

Proof. Subtracting (2.21) from (3.3) gives (5.4a) for g < |a;| < 1. By analytic 
continuation, it holds also for 1 < |a;| < q~^. The equation (5.4b) follows similarly 
from (2.7). □ 

Applying (5.4) to the matrix elements in (3.1) leads to pfaffians of the form 
pfaff {A — A* -\- B). The following lemma will be useful. 

Lemma 5.2. One has 
pfaff (oij - aji + hij) 

l<i,3<m 

^MlZl^ilf^) XI ^Sn(cT) ]^ a^(2i-l),a(2i) n ^2i-l),a(2i), 

where M is the integral part o/m/2 and {bij)Y^j^i is skew-symmetric. 
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Proof. Consider first 

1 ^ 

pfaff {ttij + bij) = ^g^(^) n'^"^(2^-l).^(20 + br{2i-l),r{2i)) 

l^''^'^'" ^ ^^^-rGSm i=l 

^ ^aF^ 1] ^S^(^)n'^^(2^-1).^(20 n ^2i-l),r(2i), 

■ sc[M] Te5r„ ies ie[M]\s 

where [M] ={1,2,..., M}. Choose a e Sm so that 

s M 

Y\_^T{2i-l),T{2i) Y\. ^-r(2j-l),T(2i) = JJ^ a(T(2i-l),CT(2i) JJ^ b^(2i-l),a{2i) , 
ieS ie[M]\S i=l i=s+l 

where s — \S\. This can be done so that, for fixed 5", r i— > cr is a bijection. 
Moreover, sgn{a) — sgn(T). Rewriting the sum in terms of a and s yields 



pfaff {(lij + h) = ^]^^f2[s) 5^ sgn(a)PJa,(2i-i),.(2i) f[ K{2i-i),a{2t). 



Replacing aij by the sign changes obtained when expanding the prod- 

uct 111=1 («a(2i-i),(7(2i) - CLa{2i),a{2i-i)) are compcnsatcd by the factor sgn((7), thus 
yielding the desired result. □ 

If ttij = Yll^=i Ck{,Xi/xjY and hij = {xi — Xj)/{xi + Xj), Lemma 5.2 gives 

(oo 
' ^ + J^Cfe Uxi/Xj)'' - {Xj/Xif\ 
Xi + Xj 

^ M ^ oo s M 

where, as before, the sum over ki, . . . ,ks equals 1 when s = 0. By (5.2), this may 
be written 

X - X ■ ^ 1 °^ * 

n ^v^xiij^ XI n^'='*5(fci,-fci,...,fc.,-M(^i'---'^'»)- 

l<i<j<m * s=0 ■ jfci,...,jts=l 1=1 

Similarly as for Theorem 3.5, we exploit the symmetry in ki to write the result as 
follows. 
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Proposition 5.3. One has 

/ oo 

pfaff ' ^ ^'^Ck ({xi/xj)'' - (xj/xi)'') 

\<^i. 'i<'m \ X'i H~ Xn _ ^ / 



l<i,j<m + Xj 

X X ^ 1 ^ 
^ n X- + ^ Ylck,Q(kr,...,ks,-k„...,-kr)ixi,...,Xm), 

l<i<j<m ' ^ s=0 ki>k2>->ks>l i=l 

where {ck)kLi is a scalar sequence such that the series converge absolutely and M 
is the integral part of m/2. 

Remark 5.4. More generally, the same proof gives 

oo \ 

?, X ^ /hi h J \ \ T T ^t, 



k,l=—oo / l<i<j<m 



^ IJ^j XI IJcfc,;, Q(fci,...,fc„zi,...,z,)(a;i,...,a;^). 

We apply Proposition 5.3 when m is even and Ck = 2{—q)^/ (1 + g^), and when 
m is odd and Ck — — 2(— — q^). Then, by Lemma 5.1, the pfafRans are of 
the form (3.1). We thus arrive at the following result. 

Theorem 5.5. If q < \xj/xi\ < q~^ for alii and j, then 



{(if^ J-j- {qxj/xi,qxi/x^ 



]JCO 



l<^k'im (-^^iM' -Q^i/Xjjoo 

"^1 * {—<iY' 

^5^2^ 5Z n 1 I ofci Q{ki,-,ks-ks,...,-k{){xi, . . .,X2m), 



s=0 ki> ->ks>l j=l ^ 

(9)^ TT {qxj/xi,qxi/x^ 



n 



(-^)~ l<^<|<2m+l (-9^.7^^ -9a^^/^j)oo 

~ 2^~^ 1^ 11 l_ofci V(fei,...,fc.-fc.,...-fci)la;i,---,2;2m+lJ- 
s=0 fei>— >A;s>l i=l 

Remark 5.6. It would be interesting to find an algebraic interpretation of The- 
orem 5.5, beyond the link to denominator formulas for queer affine superalgebras 
via modular duality. Plausibly, the work of Sergeev [SI, S2], where Schur Q- 
polynomials arise as characters and spherical functions for queer superalgebras, is 
relevant for finding such a relation. 

Specializing as before Xi = \ and replacing q by —q^ Theorem 5.5 reduces to the 
following identities. 
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Corollary 5.7. One has 

m ^ s u. 

=0 fci>-->A;s>l 



s=0 fci> ->ifcs>l j=l ^ 



□ = ^ JJ^_Q„. ,„_.„...,_,.,(l-«). 

s=0 fci>-->fcs>l i=l ^ ' 

By (5.3), this is equivalent to Corollary 4.3, and thus also to Milne's sums of 
squares formulas. 



6. A NEW FORMULA FOR 2m^ SQUARES 

As we have seen, Milne's formulas for 4m^ and 4m(m + 1) squares arise from 
the pfaffian evaluations (3.1) as all variables — > 1. In view of the results for 
triangular numbers in [Rl], one would expect that more general formulas for 4m^/d 
squares, when d \ 2m, and 4m{m + l)/d squares, when d \ 2m or d, \ 2(m + l), may 
be obtained by letting Xi tend to suitable fractional powers of q. We have made 
some preliminary investigations in this direction, but, unless one can simplify the 
arguments, the results seem very complicated. For the benefit of the interested 
reader, we state without proof a particularly simple special 2m squares 

formula which we consider to be the natural analogue of (1.5). It can be obtained 
from (3.1a) in the limit when Xi, . . . ,Xm 1, x^^i, . . . ,X2m y^- 

To state the result we need to introduce the Schur-type polynomials 

is) , ^ , deti<,,,<^(pitj-i(^0) 

jT^m \xi^ . . . , Xra) — A(a:) ' 

where p"^'' are as in §4.1. They generalize the correlation functions of §4.3, since 

Ui=l \\Pn-i\\ 

Moreover, if n — m = 2A; + £, with k a non-negative integer and s e {0, 1}, we 
have the following generahzation of (5.3) [R2, Corollary 5.11]: 

Q(Ai,...,A„)(l") 

2|m(2n+l-m) ^ -^^km ^ ^ 

lij=ll" l<i<j<m 



26 



HJALMAR ROSENGREN 



Theorem 6.1. One has 

S0,S1,'S2>0 

so+2si<m, S2<si 
2(so+si+S2)(2m+l)-i((soH-2si)2+(so+2s2)^) 

so ^0 So 

E n(-i)^'^"'^n^^'^ n 

fcl>--->fctg >1, fetQ+i>--->fcsg >1, 'i,---,'sq odd positive i=l i=l i=to+l 

fcj>--->fe^^ >1, Z'j^ even positive 

k">—>k'J^>l,l'{ even positive 
fclh+---+fc«o'-0+fcl'i+---+fc^i'n+*^l'l+---+'="2'"2=" 

X n (^1-^^')' n (^«-(^;)') n (^.-(^;)') 

l<j<jr'<to l<i<SO> l<J<Si l<*S:So, l<i<S2 

or to+l<i<j<so 

X n ((^;)' - (^D^)' n mr-ik'irr 

l<i<j<si l<*<i<S2 

X PSo+2si( — ^1, ■ ■ ■ , ~^so' ■ ■ ■ ) ~{Kl)^1 ~(^Sl)^) 

'•I 

X PSo+2«2(~^1' ■ ■ ■ ' ~^S0' ~iK)^-> ■ ■ 1 "(^sg)^' "(^52)^)' 

^^2 

where e = Q if m — Sq is even and 1 e/se, anc? where Iq, ti andt2 denote the integral 
part of So/2, {m — Sq — 2si)/2 and {m — Sq — 2s2)/2, respectively. 

As examples, we work out the cases m = 1, 2, 3. We observe that the term 
s = {sq,Si,S2) = (0,0,0) only contributes to the trivial case when n = 0, so we 
assume that n > 0. Note also that when m = 1 or 2, all Schur-type polynomials 
Pfci that appear have either A; = or / = 0, and thus equal 1. 

2 squares: When m = 1, the outer sum has only one non-trivial term, s — 
(1,0,0). We recover the two squares formula (1.1a). 

8 squares: When m = 2, we obtain an eight squares formula, which is more 
complicated than (1.1c) but still interesting to discuss. We have four non-trivial 
terms, s = (0, 1, 0), (1, 0, 0), (0, 1, 1) and (2, 0, 0), giving 

□8(n) = 16 E (-l)''^~'+^''^fc; + 16 E (-1)^^''"'^'^? 

A;iZ^=n,i^ even kili=n,h odd 

+ 64 E (-l)'i+^'i+'=i'+l'i'fc;A;'/ + 64 J] 

k[l[+k'{l'{=n k-ih+k2l2=n 
l[ and /'/ even and h odd 
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Equivalently, 



fc=i ^ fc=i ' ^ 

\fe=l ^ / fc=l ' ^ fe=l 



Computing the Lambert series using (2.10a), (2.10b) and (2.11) gives after sim- 
plification 



Applying □(— g^)^ — □(g)n(— g), which is easily verified either from (2.8) or from 
the definition, we recover Jacobi's quartic identity (2.12). 

18 .squares: When m = 2, wc have seven non-trivial terms, s = (0, 1, 0), (1, 0, 0), 
(0, 1, 0), (1, 1, 0), (2, 0, 0), (1, 1, 1) and (3, 0, 0). The second and fourth term involve 
the polynomial P^i'*(a;) = pf'\x) = x — \, while all other Schur-type factors are 
trivial. We thus obtain 



□is(n) = 25 {-if'^^HKf 

k'^l'^=n, l'^ even 

+ 2^ (-i)^(^^-^)(^i^ + ^)' 

kili=n, l\ odd 

+ 28 Y (-l)'^+^'^+'"+^'"(A;;)'(A;i)' 

l'^ and l'^ even 

+ 2« (kl + i") k[ 

ll odd and V-^ even 

+ 2« Y {-lP'^''-^^k\k? 



2 

kih+k2h=n 
h and i2 odd 



+ 2^° Y {-l)h(i^-^)+K+li'i+k'{Hi'{k'^k>i 

kxh+k[l[+k'(l'(=n 
ll odd, i J even 

+ 2^0 J] (-l)i(^i+'^+'3-3)A;2(A;2 - klf. 



kih+k2l2+k3l3=n 
hthjh odd 
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Appendix. Ono's sums of squares formulas 

In this Appendix we show that Ono's sums of squares formulas [On] are equiv- 
alent to those of Milne. Using the notation of Ono, we define by 



i=l l<i<j<m A=(ai,...,am) 



1=1 l<i<j<m \={ai,...,am) 

Moreover, we need the modular forms 



n=l 



_22k- 



2k I l-l) l^2fe| , _ /„N„2n 



n=l 



+ ^(-ir<T2,_i(n)gM, 



4A; 

n=l 

where crfc(n) = Xld|n'^''- Compared to [On], we have replaced q by —q, and used 
that = (— l)*~^|-B2i|- In generating function form, Ono's result then takes the 
following form [On, Theorem 1]. 

Theorem A. 2 (Ono). In the notation above, 

m'"^' - \ Jm-l , E ^tm-'i^l + 1) • ■■E-'i^m + 1), 
J-i-t-i A=(ai,...,am) 

□(^)4.(^+i) ^ ^ A-^{\)E-{a, + 1) ■ ■ ■ E-(a„ + 1). 

lii=l A=(ai,...,a^) 

We claim that Theorem A. 2 is equivalent to Corollary 3.6. To see this we observe 
that, in the notation of Lemma 2.1, 

^m(^l' ■ ■ ■ ' ^"i) ~ C(ki, . . . , kjn), di — '^ki + 1, 

^m(^l' ■ ■ ■ ' ^m) — C{kl, . . . , km), di — "^ki + 3. 

Thus, Theorem A. 2 can be written in Hankel determinant form as 



1 \mAm 
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Our claim would now would follow from the identities 




4(-i)'Er=^ 



Cancelling the terms involving Bernoulli numbers, we are reduced to proving that 

°° oo oo „ 2fe-l 

n=l n=l n=l ^ 

cxD cxD '-'^ n 2/;;— 1 

n=l n=l n=l ^ ' 

or, equivalently, replacing q by —q in the second identity, to the elementary iden- 
tities 

l2k-l^lm _ 2 ^ {2lfk-l^Alm ^ J- ^_^^{l-l){m-l)i2k-l^lm ^ ^^.1) 
/,m>l Z,m>l ?,m>l 

2 5^ (2/)2'=-^g2;„z _ ^ i2k-i^i^ ^ 

l,m>l l,m>l l,m>l 
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